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We present a semiclassical analysis for a dissipative quantum map with an area-nonpreserving classical limit. We 
show that in the limit h — > the trace of an arbitrary natural power of the propagator is dominated by contributions 
from periodic orbits of the corresponding classical dissipative motion. We derive trace formulae of the Gutzwillcr 
type for such quantum maps. In comparison to Tabor's formula for area-preserving maps, both classical action and 
stability prefactor are modified by the dissipation. We evaluate the traces explicitly in the case of a dissipative kicked 
top with integrable classical motion and find good agreement with numerical results. 
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I. INTRODUCTION 



A good deal of what we know about the quantum mechanics of classically chaotic systems is due to Gutzwiller's 
trace formula It expresses the fluctuation of the density of states to lowest order in Planck's constant fi by a sum 
over classical periodic orbits. Gutzwiller developed this semiclassical formalism for conservative Hamiltonian systems. 
Later on, the formalism was applied by Tabor to area-preserving maps ||. A lot is known about the quantum 
mechanical as well as the classical dynamics of dissipative chaotic systems. While the quantum energy spectra of 
conservative chaotic systems tend to have spectral statistics like the well known ensembles of random Hermitian 
matrices || , the spectrum of the quantum mechanical propagator of dissipative chaotic systems seems to obey the 
statistics of Ginibre's ensemble. This was shown numerically by analyzing a damped version of the kicked top ||. 
The eigenphases become complex and the spacing distribution P(s) of nearest eigenvalues shows a universal small-s 
behavior, P(s) cx s 3 , as soon as the dissipation is strong enough. For a classically integrable dissipative system, 
P(s) cx s. This is at variance with the result for conservative systems, with the P(s) of Wigner and Dyson, P(s) cx s* 3 
(s <C 1) with (3 = 1,2,4 depending on the time- and spin-reversal symmetry of the system. For classically integrable 
conservative systems, the spectrum is uncorrelated and P(s) — exp(— s) 0. 

Classical dissipative and chaotic dynamics typically leads to phase-space structures dominated by strange attractors, 
i.e. attractors of fractal (or multifractal) dimension, whereas for conservative systems initial distributions in phase 
space rapidly converge towards a stationary uniform one on the energy shell jg] . 

Particularly noteworthy earlier semiclassical work on dissipative processes is the quantization of Henon's map with 



i 

dissipation by Graham and Tel || ; Dittrich and Graham considered the kicked rotator with dissipation jlO| . In 1 1 1 



a 

periodic-orbit analysis of the classical spectral density for the Lorenz model, the Hausdorff dimension of the Lorenz 
attractor, and the average Lyapunov exponent was put forward. Refs. p"2fl and [ fl3| treat numerically the transition 
from quantum mechanics to classical mechanics for the kicked top with dissipation. In fj^] it was shown by direct 
numerical diagonalization of the propagator of the density matrix that the eigenstate corresponding to the lowest 
eigenvalue when written in a coherent state basis tends towards the corresponding stationary classical distribution 
function, in particular to (a smoothed out version) of the strange attractor in the case of classically chaotic motion. In 
Jk| the same analysis was carried further by solving the equivalent stochastic Schrodinger equation. This procedure 
allowed for much larger values of the angular momentum and correspondingly smaller values of Ti and showed again 
a remarkable agreement between the two stationary distributions in phase space for the quantum and classical case. 

In this paper we establish a periodic-orbit theory for a dissipative quantum map P that consists of a purely unitary 
step followed by a purely dissipative one, the dissipative one describable by a Markovian master equation. We show 
that to leading order in Ti traces of arbitrary natural powers of P are given by a Gutzwiller type formula, 

trP^£| f^=^ , N = l,2,..., (1.1) 

p.p. |tr]ltivM«-trM 

where the sum is over all periodic points of the corresponding classical map P^ in phase space. The action appearing 

in the exponent is a classical action from the dissipative steps in P^ , the matrices M^p are the monodromy matrices 
for the dissipative steps and M is the total monodromy matrix of Pzf , The inverted order of the indices at the product 



indicates that the matrices are ordered from left to right according to decreasing indices. We will discuss (1.1) in 
detail below and illustrate our reasoning for a dissipative and periodically driven top, which in the integrable case 
allows for explicit evaluation of all traces. 



By tops 15 1 we mean systems whose dynamical variables are the three components of an angular momentum, 
Jx,y,z, the squares of which have a conserved sum, J 2 = j(j + 1). In the classical limit (formally attained by letting 
the quantum number j approach infinity) the surface of the unit sphere lim^oo (J/j) 2 = 1 becomes the phase space, 
such that one confronts but a single degree of freedom. Classical nonintegrability can if desired be enforced by 
periodic driving. Quantum mechanically, one deals with a (2j + 1) dimensional Hilbert space within which one must 
diagonalize the Floquet operator to get the spectrum of (2j+ 1) quasienergies. The semiclassical limit is characterized 
by large values, integer or half integer, of the quantum number j. The convenience of periodically driven tops for 
investigations of spectral properties mostly lies in the compactness of the classical phase space and the corresponding 
finite dimensionality of the quantum mechanical Hilbert space. 

Our theory applies to physical situations where dissipation is important. One might try experimental realisations 
with several identical two-level atoms collectively interacting with a single mode of the electromagnetic field, as in 
the superradiance experiment of Haroche [p3[ . 

The paper is organized as follows: In Section [l| we introduce the map to be studied, the composition of a purely 
unitary and a purely dissipative part. Semiclassical approximations for both components of the map are reviewed. 
Section III is devoted to the derivation of a trace formula for the propagator and Section deals with the general- 
ization to higher traces. In ^ we apply the theory to an integrable dissipative kicked top and compare with numerical 
results, and in Section [v| we summarize our main conclusions. 



II. THE MAP 



While the stroboscopic description of unitary dynamics employs a wave function ip and a unitary Floquet operator 
F as ip(nT) = F n ip(0), a dissipative quantum map requires a density matrix p{nT) — P n p(0) with a nonunitary 
propagator P. For the sake of convenience we imagine dissipative and purely unitary motion to take place separately 
and in sequence during each period such that the map takes the form 

p((n + 1)T) = DFp(nT)F ] = Pp(nT) ; (2.1) 

in contrast to the unitary part of this map (which factorizes in F and F') the dissipative part D cannot be split into 
two factors. Our goal is to determine spectral properties of the total propagator P, a linear, non-unitary operator, 
by semiclassical means. As a preparation we briefly review what is known about semiclassical approximations of the 
two components of the map, starting with the unitary part. 



A. Semiclassics for the unitary evolution 

The Floquet operator F may by itself describe a motion with either an integrable or chaotic classical limit. A simple 
integrable example is a pure rotation by an angle (3 about, say, the y-axis with F = exp(— \(3J y ). Most illustrations 
of general results will in fact be given for that case. If desired, a chaotic top can be treated by including a torsion 
about the z-axis as F — exp(— i 2 J 2 ) exp(— i(3J y ). The corresponding classical motion proceeds in the spherical 
phase space defined by the conservation law J 2 = j(j + 1) and may be described by a pair of canonical variables; a 
convenient choice is attained by specifying the orientation of the classical unit vector lim^oo j by a polar angle O 
with respect to the z-axis and an azimutal angle <fr reckoned against the x-axis in the xy-plane and then taking cos 9 
as the momentum and <fi as the conjugate coordinate. 

A semiclassical approximation for F can be derived by representing F in the J z basis (J z |m) = m\m)). The matrix 
element F nm = (n\F\m) is the probability amplitude to get the z-projection of the angular momentum mapped from 
the initial value m to the final one n. If F describes a pure rotation that matrix element becomes Wigner's d-function 
whose asymptotic behavior for large J = j + \ is well known |l6j| . The discrete quantum numbers m and n can be 
replaced in the limit J — ► oo by the continuous classical initial and final momenta 

p = raj J = cos 0j, and v = nj J = cos 9 / , (2-2) 

respectively. The semiclassical approximation has the familiar van Vleck form in which a classical action S(u, p) enters 
as a phase, and its second mixed derivative 

C(u, p) = { - 7 =tL j\d^S{v,p)\ (2.3) 
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as a preexponential factor. All classical orbits connecting the initial and the final momentum contribute a term of 
the structure G{v, p)e lJs ^' l ' l \ Explicit expressions for the function S{v,p) can be found in Esl, both for integrable 
cases such as pure rotations and chaotic ones involving rotations and torsions. For our present prime example, the 
linear rotation by (3 about the j/-axis, there are typically two classical trajectories differing by the signs of azimuth 
and action; only if (3 > 9i + 9f or (3 < \9i — 9f\, no classical trajectories connecting the initial with the final momentum 
exist; if one of the foregoing inequalities is turned into an equality we have <pi = <pf =0 and only one trajectory 
contributes (see the discussion in |l8|]); we shall not consider these latter exceptions here but only treat the case of 
two contributing trajectories for which the amplitude under discussion reads 

F nm ~ C{v,ti)e^ Js ^-V . (2.4) 

a=±l 

The action S(v, p) (or — S if <j> < 0) also serves as the generating function of the classical map 

(M,<M^(y,4>f) (2-5) 

in the sense 



(2.6) 



here, the upper (lower) sign refers to <f>ij > (< 0); note that for the rotation about the y-axis the sign of the initial 
and final azimuth is the sa me. 

The generating property (2J3) and the general van Vleck form (2^4) will allow us to develop semiclassical formulas 
largely independent of the particular dynamics. In particular, the specific form of S(v, p) will not be needed below. 



B. Semiclassics for the dissipative part of the map 

To introduce the dissipative part D of our map ( |2.l| ) we momentarily allow for a continuous time and employ the 
master equation of spin damping well known from the theory of superradiance 



p{t) = K ([J_,p(t)J+] + [J-p(t),J+]) = A P (t) , (2.7) 

where the damping constant k sets the time scale on which the angular momentum in question approaches, while 
keeping its square J 2 = j(j + 1) conserved, the stationary state of minimal ^-component, J z = —j; the operators 
J± = J x ± i J y are the familiar raising and lowering operators. 

Let us first look at the classical limit. The classical equations of motion can be determined by extracting equations 
for the expectation values of J z and J±, parameterizing them as (J z ) = JcosG and (J±) = Jsin^e* 1 ^ and factorizing 
all operator products (e.g. ( J+ J_) — > (J+)(J_)). One then finds that the angular momentum behaves classically like 
an overdamped pendulum |lq ], 

= 0, e = 2jKsin6. (2.8) 

The classical damping constant is thus revealed as 2 Jk. For future convenience, we give the solution of ( p.S| ) in terms 
of the dimcnsionless time s = 2Jnt and phase-space variables p(s) — cos0(s) and <fi: 

1 (i-Ms))(i + /i(Q)) , , , „<. 

s= o ln 7i tnwii i — nvT' 4>{s) = 4>{0) ■ (2.9) 

2 ( 1 -A*(0))(l+M(s))) 

For the sake of definiteness we imagine the damping to be effective during the full duration T of a period, thus 
confining the previous unitary motion to an instantaneous kick. The dissipative part of the map is then obtained by 
exponentiating the continuous-time generator A of the above master equation as D = cxp(AT). 

The semiclassical approximation for the dissipative propagator D = exp(AT) looks much like the familiar van Vleck 
form for unitary propagators. Details of the theory which amounts to solving a Hamilton-Jacobi equation for the 
classical generating function of the dissipative motion can be found in |2~ijj . Here we limit ourselves to exposing the 
idea and summarizing the main results. 



Looking at the master equation (ETtT) in the J z basis we reveal an important conservation law upon expressing 

^sa and the skewness k = ^ 



the density matrix element (mi\p(sj\rri2} as a function of the mean m — mi + m2 and the skewness k = mi 2 m2 of 



3 



the quantum numbers mi,m 2 , i.e. p m (k,s) = (mi\p(s)\ni2)- The operators J± give rise to the functions d„ 



3 \i + 1) — m(m — 1) known from the angular momentum algebra ( J± \m) = ^/d^m\rn± 1)) and f m (k) — y/ d m +kd m -k- 
With these abbreviations our master equation reads 



j dp m (k, s) _ j m (fyp m+1 (k,s) - (d m - k 2 )p m (k, s) . 
as 



(2.10) 



Obviously, the skewness k enters only as a parameter and is conserved. The general solution of (2.10) can be written 
with the help of a "dissipative propagator" D mn (k, s) as 



p m {k, s) = 22 D mn (k, s)p n (k, 0) 



(2.11) 



For large J, the variables p — m/J and v = n/J prove useful again. A third one is introduced for k, T) = k/J, and 
we write consequently D mn {k, s) — J~ 1 D(p, v; ij; s). 
The clue towards D is a WKB type ansatz, 



D(p, v; rj; s) = B(p, v; r]\ s) exp(JR{p, ij; s)) , 



(2.12) 



and to expand systematically in 1/J. Note that even though the functions B and R are assumed smooth, the WKB 
form still allows for a delta peak as J — > oo. To lowest order in 1/J, (2.11) leads to an equation for R which looks like 
a Hamilton-Jacobi equation for a one dimensional Hamiltonian system with the somewhat unconventional Hamilton 
function H(p,p) = (1 — p 2 ) (1 — exp(p)), 



OR dR. n 

-ds- + H ^dp- ) = °- 



(2.13) 



In this Hamiltonian description of the semiclassical approximation of our damped spin the v ariabl e p appears as a 
coordinate with a canonically associated momentum P = f^- The Hamilto-Jacobi equation ( 2.1 3| ) is equivalent to 
two Hamilton equations of motion for p and p. The solutions involve two integration constants, one of which we 
choose as the initial coordinate v\ as the second will serve the conserved "energy" E — H(p,p) which for the sake of 
convenience we represent by the parameter a = y/1 — E. With this choice, we get the Hamiltonian trajectories 



1 (v + a)(p-a) 
— m — - 

2a (v — a)(p + a) 



(2.14) 



Comparison with the classical dynamics ( |2.9| ) shows that p and v are connected by a classical trajectory iff a = 1. If 
one thinks of a Hamiltonian trajectory as specified by th e initial and final coordinates v and p with a given time span 
s, as we shall frequently do below, one should read (2.14) as determinin g the parameter a as a function a = a(p, v, s). 
The action R can be obtained by integrating pdp along the trajectory ( 2.14 ) and reads J2l| 



R(p, is; ry; s) = -(£(1, v - rj) - £(1, p — rj) + £(1, v + rf) - £(l,p + rj)) 

-Z{ ai is)+Z(a,p) + s(a 2 -1 + f] 2 ), 
£(x, y) = (x + y) hx(x + y) - (x - y) ln(x - y) . 



(2.15) 



Several features of the action R will be of importance in the sequel: 



• For rj — we get a — 1 from d^R = 0, thus the classical equation of motion (2.9); then d v R = holds as well. 

• For r\ = and p, v connected by the classical trajectory (i.e. a = 1), R is strictly zero. This actually holds 
beyond the semiclassical approximation, as it can be traced back to conservation of probability by the master 
equation for k = po[ . 

• R is an even function of rj and has in fact always a maximum at 77 = 0. 

An equation like ( 2.13| ) has already been found for a dissipative problem in Jl9j , but the preexponential factor which 
is crucial to our problem was not treated there. That prefactor is obtained in next-to-leading order in 1/J as [|2~i| 



B(p,is;r);s) = 



\J d^d v R{p, v\ rp s) ; 



(2.16) 
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it differs from the familiar form */ d^d^Sjy, p) in the van Vleck formula for unitary propagators by the additional factor 



| Q . Since ^Jd^d^S^, p) arises from the propagation of a wave-function we might expect the square of this factor 
when propagating a density matrix with r\ = (i.e. probabilities). And indeed, both square roots in ( 2.1 6| ) are the same 



in the semiclassical limit and combine to the classical Jacobi determinant: For J^oowe can restrict our attention to 
values of p, and v close to the classical trajectory, i.e. a = 1, as the propagator for other values is exponentially small. 

We call this trajectory pd{v) where the time is fixed. The first square root just gives J d>i %il ^ n ^ ne °t ner hand, 

since R(p, v, 0; s) has a maximum on the classical trajectory, it can be written as v; 0; s) ~ ~§(^ — / i d 1 (/^)) 2 

close to the maximum with some parameter a > 0. This means d^d u R{p, v] 0; s) — a -g* , and for J — > 00 ( 2.16| ) 
therefore leads to the classical propagator of probabilities, 

d V~4 



D c i(^i, v, 0;s) 



dp 

As a conclusion, we may write the preexponential factor for r\ — as 




(2.17) 

if it is evaluated on the classical trajectory. This is what we will need for the semiclassical evaluation of the traces of 
P (see the next section). 

A word is in order here about the limits of validity of the present semiclassical treatment of dissipation. It disregards 
the existence of discrete energy levels and becomes unacceptable when this discreteness is important. These are the 
cases when 

• the initial state n and the final state m coincide or are separated by just a few intermediate levels: n — m ~ 1; 

• m or ti (or both) coincide with the highest or lowest energy state of the system | ± j) or are separated from 
them by just a few levels, i.e. when at least one of the following relations hold ±j ± n ~ 1, ±j ± m *~ 1. 

The first limitation implicitly puts a lower bound on the time, s 3> 0{J~ l ), because for small propagation times 
it is only the elements D mn with n — m zero or of order unity which are noticeably different from zero. The failure 
of the continuum approximation is self evident in this case: the probability literally had no time to get away from 
the initial level to form a smooth distribution. The second limitation means that the initial and final coordinates 
should not be too close to the "poles" ±j. The inapplicability of the naive semiclassical approximation to describe 
the lowest energy state of a quantum system is well-known; in systems with the spectrum bounded from above this is 
also true about the highest-energy state. Again, there is an implicit limitation on the time which should be not too 
large (small compared with J) so that the probability is not concentrated close to the ground state. 

In the present work the deficiencies of the semiclassical approximation just mentioned seem to be of no importance, 
at least for the lowest traces of P. For higher traces we see deviations to next order in 1/ J which may well be due to 
our neglecting the discreteness of the spectrum. A better approximation for the dissipative propagator exists which 
allows for uniform evaluation, at least for not very large times EJ. 

Since both F and D have a classical counterpart, so does the combined map P. We shall denote it as P c i and 
propose to show that its periodic points determine the spectral properties of P. 

III. SEMICLASSICAL EVALUATION OF THE FIRST TRACE OF P 

We first write out an exact formal expression for the quantum map P in terms of the discrete quan- 
tum numbers. Immediately after the unitary motion, p has the elements j o„(fc,0+) = (ni\F p(0)F^\n2) = 
Yj > i u i 2 ^ 3 Fn 1 i 1 {F^)i2n 2 {h\p{m2) = Ei,r F n+k,l+ r FZ_ k<l _ rP i(r,0) with / = ii±ii and r = We insert this 

in the general solution ( |2.1l| ) and get p m (k, r) = Yin It D mn{k, r)F n+ k,i+ r F*_ kl _ r pi{r, 0) = J2i r p mk-irPi{r, 0) where 



we have introduced the dimensionless parameter r = 2JkT. The propagator P is read off as 

P-mk-lr = ^ Dmn(k,T)F n+ kA+rF*_ k j_ r , (3.1) 
n 

and yields the first trace 



5 



trP = /] Pmk;mk = D mn {k,T)F n+k ^ m+k F*_ km _ k . (3.2) 

m.fc m.n,k 

We now replace sums by integrals via Poisson summation, employ the semiclassical propagators derived in the previous 
section, and shall eventually integrate by the saddle-point method. The integrals to be done read 

OO r. 

trP = J 2 / 2dfxdvdr]exp[J^(fj,, v-rj)} ^ 

r,l.t — — oo <J\ ,&2 

■B(fj,, v, rj)C(v + r), fi + rj)C* (v - rj, fi - rj) . (3.3) 

They involve the now complex action 

$(/j, v; rj) = R(p, v; rf) + i27r(r(p + rj) + l(v + rj) + 2trj) 
+\{<7iS{v + ?], f-i + rj) + o 2 S(y — r},fi — r/)} 

= $ (M,^)+i$iG",z';»7) (3-4) 

whose real part is related to dissipation while the imaginary part represents unitary dynamics (We apologize to the 
reader for frivolously calling action what is usually called — ix action). The integral contains factors 2 in front of 
the integer t and before d/i which arise because m, n and k can be simultaneously half-integer. For brevity, we have 
suppressed the parameters /3 and t. 

Three complex saddle-point equations (SPE) have to be fullfilled: (9 M $ — 0, 9„ < I> — and = 0. We look for real 
solutions (complexity would indeed look unphysical, as m — jij G Z etc.). Formally we can not exclude the existence 
of further complex solutions, but as long as classical solutions exist we expect them to dominate over non-classical 
ones. This is certainly true in the case of the superradiance dissipation where R = on the classical trajectory, 
whereas complex solutions would lead to exponential supression from the term i<f>i. But even more generally we 
expect classical solutions to dominate, since they are known to dominate in non-dissipative quantum mechanics and 
dissipation favors classical behavior even more. 

We thus need separately = and d^&i = 0, where <&o and i$i are the real and imaginary parts of $ (and 

correspondingly for the derivatives d v and d v ). These are six equations for the three variables /i, v, rj, but the additional 
three integers r, I, t from the Poisson summation do allow for a solution. 

The first equation, d v &o — 0, immediately gives i] — with the important consequence that only the propagation 
of probabilities contributes to the trace in the semiclassical limit. Inserting 77 = in d v $>i = 0, we are led to 
(a 1 -<T 2 )(d u +d l ,)S(v,tM)+2n(r + l + 2t) = 0. 

Consider first the case u\ = —<J2- In view of the generating properties of S{v, ^),the foregoing equation yields 
(j)f(v 7 [i) — (pi(u,y) + n(r + I + 2t). From d^i = = 2nr + (tri + a 2 )d fl S(u 7 n) = we conclude r — 0; and from 
d v <&\ — correspondingly / = 0. We thus obtain the unique classical trajectory for the unitary evolution which 
connects the initial and final momenta for the given time span r, uniqueness due to the additional requirement that 
initial and final azimuths be the same up to a multiple of 27r! The second of the three saddle-point equations for 
the real part, <9„$o = 0, immediately leads to a = 1 (since 77 — 0), i.e. precicely the condition distinguishing the 
trajectories of the classical overdamped pendulum within the larger family of Hamiltonian tracjectories describing 
the semiclassical limit of the dissipative motion. With that condition met, the last equation = c'^'l'o = — du&o also 
holds without further restriction. 

To summarize, the saddle-point equations determine the classical trajectory which leads, in the unitary part of the 
map, from /1 to v with constant <j) and, for the dissipative part, back from v to /1, again under constant (f>. The saddle 
points of the above integrand are the fixed points of the classical map. 



Consider now o~\ = a 2 = o. The imaginary part of the rj variation, d v ^>\ = 0, does not restrict /z, r\ but just gives 
r + / + 2t = 0. New information now comes from 9„<&i = and 9 M $2 = and amounts to <fif = ttI = — 7r(r + 2t) and 
(pi = —nr. So we have again (pf — up to an integer multiple of 2-7r, but additionally (pi = or 7r. The dissipative part 
remains unchanged, such that for the combined map we get the particular class of periodic orbits that take place at 
(p = or 7r. For our integrable dissipative kicked top, such orbits are only possible if the damping constant r exceeds 
a certain critical value r c . Moreover, we shall see that these orbits lie exact ly on the boundary between classically 
allow ed and classically forbidden motion concerning the rotation alone (cf. sec II A). The semiclassical approximation 
(2.4) then ceases to be valid and should be replaced by a uniform approximation bridging the classically allowed and 
forbidden domains |l8| ]. This program is beyond the scope of the present work; we will restrict ourselves to regimes 
where this particular kind of orbits does not exist. 

Note that the combination o\ = —<j 2 would arise automatically, if the van Vleck propagator contained only one 
classical path. The opposite sign results from the time reversed propagator . It means the interference of a given 
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orbit with itself. 



The fact that only the interference of an orbit with itself contributes can be seen more generally in the case where 
the unitary part contains several classical paths with actions Si{v, /x). One encounters then a term iJ(Si(v + 77, /i + 
77) — Sk{v — T},fi — 77)) and a double sum over paths I, k. However, the SPEs lead to an equality of initial and final 
angles (modulo 2tt) for the two paths I and k. Since to a given initial phase-space point there is only one classical 
path, only I = k can contribute. Even the case of I and k designating two branches of a separatrix in phase space can 
in general be excluded, due to the equality also of the final phase-space point. The only exception are degenerate 
paths by which we mean that two originally separate paths conincide when a system parameter is suitably changed. 
Exactly this extraordinary situation occurs in the above example for the kicked top: For r < r c there are always 
two symmetry related classical paths for the unitary part (parametrized by a = ±), one starting at an azimuth </>,, 
the other at an azimuth —<pi. When r is increased both conincide for r = t c at (pi = and remain there for all 
values r > r c . But at the same time classical dynamics breaks down and the van Vleck propagator ceases to be valid, 
anyway. 

To exploit the SPA, we have to expand $ to second order in the deviation x T = (S/j,, 8v 7 Sr/) from the fixed point 
(/iO; u o] 0). We obtain a quadratic form with a complex symmetric 3x3 matrix Q3X3 that contains the second partial 
derivatives d Xi d Xj $>((j,, v; rj) L „ . in an obvious way, 

v; 77) ~ $(/i , i>o; 0) - ix T Q 3x3 x . (3.5) 

We have chosen to pull out a minus sign from the partial derivatives. The fact that Q3X3 is complex and therefore 
non-Hcrmitian makes the SPA non-trivial. In view of the repeated application of the saddle-point method for the 
integration of a complex function of several variables, we devote appendix ^ to a discussion of the general formula 
derived in |2l|J27|]. 



Let us come back to our saddle-point integral (3^). According to equation ( |Al| ) we have to evaluate detQ3 X 3. It 
is easy to see that the 2x2 (fi,v) block and the (77,77) block of the imaginary part of Q3X3 are zero. In addition, 
dpdri'&o = = d„dri&Q, since 4>o has a maximum at 77 — for all /z and u, i.e. it must have the structure $0 — 
fo{n, v)r^ for small 77 with some function f . The second derivatives of the real part can be simplified on the classical 
periodic orbit: Since d M $ (n,ir,a) = for all v, we find |-9 M $ (/j, v\ 0)^=^) = ^$0^77 + d^d^o = 0, 

i.e. d^du^o = -dfoo^-. Similarly, starting from <%$ \ ll=lld (v) = °> onc nnds d l®° = d l®v (^r) • We introduce 
the functions h{v, (i) = 2(d u + d^)S(v, /1) = d v §i and 7 = to write Q3X3 in the form 



Q 3x3 = - I -7^*0 7 2 ^1>o id v h I . (3.6) 



7 2 af<i 

id u h \d v h <9^$o 



It has the determinant 

detQ 3x3 = -dl$ {d v h + jd»h) 2 . (3.7) 

Since at the saddle point d^&o < (maximum of the propagator on the classical trajectory!), we have detQ3 X 3 > 
and by the same reason the first minor D\ = — 9^$o > 0. The second upper left minor gives zero, but we can fix this 
by introducing a small positive e in the first element. This leads to D-j = —ej 2 d^&o > 0. Thus, all minors are real 
and positive so that argdetQ3 X 3 is now uniquely fixed to zero (cf. ( |A5| ) ) . 



y/dfj,d v R{iJ.,iy; 0) = 



"' " ''" ' - — ^^fc^$o = «/— 9^$o- From the action in the exponent only the real part R — R(n,v;0) remains. 



dfi 

We arrive at the saddle-point approximation for trP 



p.o. 



dfj,d„S(v, fi) 



e JR . (3.8) 



The sum is over periodic orbits of the classical map and it is understood that all quantities are evaluated on these. 
The foregoing trace formula can be simplified further if we go over to phase-space coordinates /1 and <j> = </5j and exploit 
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the generating properties of S(is,fi). Let us call i/ r (n,<j)) the momentum component of the classical trajectory for the 
unitary motion (the index r stands for rotation), if the initial values fi and <f> are given. The corresponding angular 
coordinate will be denoted by <p r (fj,,ip); it is related to 4>f{v,^) by (f>f(v,n) = <p r (/j,,<j)i(u, jj,)). The corresponding 
quantities for the dissipative motion are \idiy, 4>) = f-div) and <^<j(f, 4>) — 4>. The generating properties of S(v,n) 

reads 



trP = V : — , „ . (3.9) 



90 di/ c? u 



Finally, the preexponential factor can be written in terms of traces of classical monodromy matrices. The monodromy 
matrices for the dissipative and the unitary motion are 

•U,= | • „ t . , I Mr=\ (S<* I • 




respectively. Their product M = MdM r is the monodromy matrix for the combined map. Recall that the dissipation 
conserves 0, such that ^7 = U- We therefore arrive at 

EpJR 
p o \trM d - trM\ 1 > 

This is the first central result. It generalizes Tabor's formula for classically area-preserving maps to an area- 
nonpreserving map and shows that even in the case of dissipative quantum maps, information about the spectrum 
is encoded in classical periodic orbits. All quantities must be evaluated on the periodic orbits and with 77 = 0. We 
remark that the formula holds for both chaotic or integrable maps, as long as the periodic orbits are sufficiently well 
separated in phase space such that the SPA is valid. 

For comparing (3.11) with Tabor's result, tr_F = ^2 p ^-trMl 1 / 2 ' one snou ^ remember that we consider the prop- 



agator of the density matrix, but Tabor the propagator of the wave functi on | a|. In the limit of zero dissipation, 
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we should get tr_P = |tr_F| . That limit can unfortunately not be taken in (3.11), since our semiclassical dissipative 



propagator is only valid for r > 1/J. However, |trF| 2 for r — ► would definitely lead to a double sum over periodic 



orbits. Of the double sum only a single sum remains in (3.11); all cross terms between different orbits are killed by 



dissipation (this is due to 77 = 0!). An important quantum-mechanical effect of the dissipation is indeed the very 
rapid destruction of interferences on time scales much shorter than the classical time scales fl5| . Decoherence leads 
in the trace automatically to the "diagonal approximation" which suppresses interferences of different orbits! 
Tabor's preexponential factor is reproduced in the limit r — > 0, since Md becomes the unit matrix, and thus trMd = 2 
and trM = tr(MdM r ) = trM r . It is, however, raised to the power 1 instead of 1/2, since we propagate a density 
matrix and not a wave function. The action itself is purely real; the imaginary parts from the unitary motion cancel 
each other at 77 = 0. This property holds for all dissipative processes describable by a Markovian master equation 
that leads to a dissipative propagator with a single maximum (as a function of rf) at n = 0. As outlined in the 
section on the dissipative part, the saddle-point action is zero in our case of spin damping. Then all physics is in the 
preexponential factor, which describes the stability of each orbit! If R ^ 0, only a single periodic orbit may dominate 
the trace in the limit of J — > 00, namely the one with the maximum value of R. 



Before we calculate higher traces, let us comment about the generality of our result. In deriving (3.11) we made 
use of the general van Vleck forms of the propagators and the generating properties of the actions appearing in them. 
Within the class of dissipative processes considered the only ingredients part icular to the specific problem at hand were 



the fact that the dissipation conserves <fi and the sum over two paths in (2.4). However, since the dissipative pro cesse s 
considered impliy only interference of each orbit with itself, the latter feature does not restrict the generality of (3.11). 
The formula should hold independently of the number of classical paths that contribute to the unitary motion. The 
vanishing of the Maslov phase relies only on properties of the dissipative part and will always hold if the dissipative 
propagator has a maximum on the classical trajectory. 
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IV. HIGHER TRACES 



We will now show how the semiclassical trace formula ( 3.11 ) generalizes to higher traces trP N . For our 2j + 1 
dimensional Hilbert space, the knowledge of (2j + l) 2 traces suffices in principle to reconstruct the entire spectrum. 
As a first step we look at the second trace. 



A. The second trace 



In principle, trP 2 could be calculated just as trP by writing down the exact discrete sums as in ( |3.2|) , transforming 
them to integrals, and doing the integrals by the SPA. The quadratic form in the exponent would then give rise to 
a 6 x 6 matrix, and in general, for the JV-th trace, to a 3N x 3iV matrix. It turns out that the problem simplifies 
considerably if one integrates over v once for all in the propagator itself, so that the latter will only depend on actual 
initial and final coordinates. The matrix in the quadratic form arising from the SPA is then only a 2N x 2N matrix. 
Let us therefore write an arbitrary matrix element of P as 



mk;m' k' 



C*(u - r),n' - 7j')exp {JG(n,r}',n',rf; v)) , 



l — — oc 



where again m = fiJ etc. and C(/i, v) is defined in (2.3). The complex action G reads 

G(fi, T); //, r)'; v) = R(fj,, v\ 77, r) + \a x S{v + 77, fj,' + 77') + ia 2 S(v - rj, // - T}') + \2irlv . 
The saddle-point equation 

d v G = = d v R + i(aid v S(v + rj, fj,' + rj') + a 2 d v S{v - v , // - 77') + 2nl) 



(4.1) 



(4.2) 



(4.3) 



will in general have complex solutions v — 77(/i, 77; //, 77'), as long as 77 and 77' are not zero and \x and p! are not connected 
by a classical trajectory. It may even be possible that several solutions or no solution at all exist. However, note 
that the saddle-point equation (4.3) simplifies considerably and gives a physically meaningful fixed point if o\ = —02 
and 77 = 7/ = I = (and that is the situation we encounter for our dissipative integrable kicked top): We then have 
d u R{p, v; 77, t) =0 which is equivalent to the classical dissipative equation of motion T> = /i^ 1 ^) and a = 1. 
When the v integral is done by a SPA the second derivative <9 2 G|„ = 77 = c£G comes into play; it can be combined with 
the other preexponential factors of P. One needs a relation bet ween second derivatives of G which can be obtained 
by differentiating with respect to fi in the saddle-point equation (4.3), accounting for the \x dependence of V, 



€g = -d^-^ 



'd^R^i, v\ T], t) 



1 



(4.4) 



With the abbreviation ^(//, 7]; \jl 77') = G(/i, 77; fi f , 77'; ^(/i, 77; 77')) for the action at the saddle point we find 



E EE 

1= — OO <7l ,(72 V 



' dv I dv 
<9/i y dp 

■C(V + 77, n' + i/)C* {V - 77, // - 77') exp ( J^(p, rj; //, 77')) 

The sum over 17 picks up all relevant saddles. To avoid possible confusion we should add that of the t wo d erivatives 
^ only the one originating from the semiclassical dissipative propagator is meant at constant a, cf. (2.16), as long 



(4.5) 



as fi and // are not connected by the classical dissipative trajectory. 

We are now ready to apply the map twice. Let us call R l and D l the i-th unitary and dissipative step of the iterated 
map. The coordinates (/i, 77) are transformed in the following way: 



n 1 



D 2 



(pwni) {vi,m) — > (^2,772) — ► (v2,m) — ► 0*3, m) = (f*i,m)- 

The last equation is the periodicity condition for trP 2 . The total action ^2 that appears in the expression for trP 2 is 
basically the sum of two actions <I> (cf. ( |3.4| )), one from each repetition of the map. All variables carry an additional 
index i = 1, 2 counting the iteration of the map (in the case of the a it is the first index), 







* 2 = R(jii,V2',rii) + R(fi 2 ,Vi;r] 2 ) + \a 2 iS(v 2 + 771,7*2 + m) 

+ia 22 S(V 2 -771,7/2 -772) + i<TuS(Vi +r) 2 , fii +i<7 12 S(Vi -772, Mi — »7i) 

+i27r((ri(/Ji + 771) + h(V 2 + 771) + 2*1771 + r 2 (fi 2 + rj 2 ) + h(pi + 772) + 2t 2 ri 2 ) . 

The last line comes from the Poisson summation over 7*i, Zj, and U. The /z's and r/'s are the independent integration 
variables whose values will now be determined by saddle-point equations. 

We assume again that the solution of the SPE's is real and split them into their real and imaginary parts. From the 
real part of the SPE's with respect to 77, d m ^ 2 = = <9 1)2 , J , 2 , we get 771 = rj 2 = 0. The real part of d fll i & 2 = 9 M2 *2 = 
leads to the classical dissipative equations of motion at constant <fi. The imaginary parts of d m i S 2 = = d ri2 1 $> 2 give 

(<t 2 i - (T 22 )(j)f(V 2 , fi 2 ) = {an - cri 2 )(j)i(Vi, fj,i) + (n + h + 2ti)2n (4.6) 
((T21 - a 22 )(j>i(V 2l ju 2 ) = (011 ~ ci2)<A/ (^liMi) + (f*2 + h + 2t 2 )2n. (4.7) 

These are conditions on the initial and final azimuthal angles for the two classical maps corresponding to the two 
unitary parts. Let us look at them more closely. Suppose that in any of these equations an — a^. We would 
get fa = or <j)f = 0, and we have excluded such cases from our present analysis. If on the other hand we had 
(<T2i — 022) = —(on — 012) 7^ 0, this would imply <f)f(Vi, fii) = — 4>i(p 2 , ijl 2 ) (modulo 71"), which is not possible since 
during dissipation (j) is conserved. We therefore conclude that only the two combinations an = —aa = a = ±1 
contribute. Then the SPE's define a period-2 point of the combined map. 

If we expand ^ 2 around a periodic point into a quadratic form, a 4 x 4 matrix Qix4 of second derivatives appears 
which consists of 2 x 2 blocks (74,17), (77,7*), and (77,77) where (/j, 77) contains the mixed derivatives —d^ li d rij ^ 2 (and 
correspondingly for the other blocks): 



74x4 — 



(/j, 77) 
(77, n) (77,77) 



(4.8) 



The {Hifi) block vanishes. To see this, observe first that its imaginary part gives zero as an — —ai 2 . For the real 
part we have by construction d-y 2 ^> 2 = = dj7 2 R(/j,i,V 2 ;0) for all fii. But then also d^ 11 R{^i,v 2 ]0) = holds for 
all Hi (cf. the general properties of R\). This immediately gives c^c^vE^ = 0. A corresponding argument holds for 
d^ 2 dfj, 2 ^ 2 . Furthermore, as 77 2 = (j,^ (fii) on the classical trajectory, 7x1 does not even appear in di_ ll R(ni,T' 2 : i 0) and 
we thus obtain d ll2 d fll R([j,i,T' 2 :0) = 0. 

Given the structure of the matrix Q4 X 4, the (77,77) block is irrelevant, since p8[ ] 

dctQ 4x4 = (det(7M7)) 2 . (4.9) 

In order to calculate the mixed derivatives in the (/j, 77) block, one needs partial derivatives of the F's with respect to 
the 77's, as we have for example 

d m d^ 2 = d U2 d tll R^ + 2iad fll d Ul S(v ut ii); (4.10) 
0772 

they are obtained by totally differentiating dn i ^ = with respect to the 77's and then setting 771 = 772 = 0; the one 
needed for the above example reads 



dv2 1 _„. dj7 2 df, 2 S {v 2 ,[i 2 ) 



2i* fl2 £ - , ■ ( 4 - n ) 



Upon transforming the remaining second mixed derivative of R in ( f .10] ) into second pure derivatives as in the case 
of trP we can eliminate the derivatives of R entirely. We find for the above example 



d m d^ 2 = 2ia\^^d T ; 2 d^ 2 S{v 2 ,H 2 ) + d^ 1 d Ul S{vi,Hi)'j . (4.12) 

In a similar way the remaining three partial derivatives of ^ 2 can be found. If we abbreviate Si — <Si(F;, fii) (i = 1, 2, 
V 3 = Vi) and Ri = i^(/ij,I7 i+ i; 0)), det Q ix4 , is given by 

detQ 4 x4 = 4 2 [ ^ ((dZ 2 S 2 )(d? l2 S 2 ) - (ch 2 d, 2 S 2 f) + p- m^idlSi) - ((h^Si) 2 ) 

+ ££ i^ s ^ s ^ - (^^SiH^A^)) 

+ (dlSi)(dl 2 S 2 ) - (dp^SiXdp^Sz) ] 2 (4.13) 
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and entails 



2 ^ lidu^S^d^S^^ 
tiP z =4 



^ |detQ 4 Jl/2 



e J(R 1 +R 2 ) _ (4-14) 



^4x4| 



The sum is over all periodic points of the map P 2 . The various prefactors from the dissipative parts have canceled, 
up to two factors g^* 1 — d^ d +1 = V ^dr"'^ . Again, for each periodic point the overall sign of the summand should 
be determined. It turns out to be positive (for a proof we refer to the discussion in the case of the iVth trace below). 
The preexponential factor can again be written in terms of elements of the classical monodromy matrices, as can 
be verified in a straightforward but lengthy calculation. We introduce monodromy matrices and Mr f° r the 
dissipative and unitary evolution in the i-th step, respectively, and 

M(') = MfpM® (4.15) 

for the entire i-th step. The total monodromy matrix M is given by the product M = M"W . . . AfW = rj| =Ar M«. 
With this convention, the second trace reads 

tvP 2 = V - . (4.16) 



B. The 7V-th trace 



Our results for the first and the second trace suggest a generalization for the 2V-th trace, 

j y r 

tvP N = V i , 6 ' , (4.17) 



p.p. |trn!= N M f- trM 

where the sum is now over periodic points of P£? . The rest of this section is devoted to the proof of the formula. 
Starting point is the intermediate semiclassical form ( |4.5| ) of the propagator P. We have to consider a sequence of 
intermediate coordinates of the form 



(4.18) 



Vl\ R 1 fV>2\ R^ (V2\ f^i\ R N f \ Df^ C^N+A _ f Hi 

mj Vw \m) Vw vw xm+i) Kvn+J \m 

All variables are periodically continued in their indices (vn+i = ^i, etc.). The total action is given by 

N N 

= '^R^^Vi^rii) + (<rnS(Vi + r/i+i^M. + Vi) + <^i2S{v l - rj l+1 ,^i - rji)) . (4-19) 

i=l i=l 

We have omitted for clarity the term coming from the Poisson summation which is easily seen to lead to integer multi- 
ples of 2n in the equality of initial and final angles that we shall find presently. For convenience we introduce a vector 
notation for the sets of variables, fi = . . . , /xjv, and similarly for V and r\. The functions Vi = rji+i', Hi, rji) 

satisfy 

d T7 ^ N {n,V,r ) ) = Q. (4.20) 

Everything goes through for the saddle-point equations as for the first and second trace. Thus, an — —oa = a = ±1, 
and the SPE's pick up the periodic points of P*[ . The 2N x 2N matrix Q2NX2N has the same block structure as in 
the case N = 2, that is d^d^^ = for all Therefore its determinant depends only on the N x N (/i,?y) block 
which we will call E>n x n in the following [ p8| : 

det Q 2Nx2N = (-l) N (det B NxN ) 2 . (4.21) 

The partial derivatives in SjvxJV are determined in the same way as for N = 2 and contain only derivatives of the 
imaginary part of \P. We find 
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d^d^N = 2ia(d ! ; l _ 1 d n _ 1 Si-iSi- 



i,k + 



Si,k 



<9 Mi Si- Si+ik 

dfii+i 



(4.22) 



where Si k is the Kronecker-delta. 



After the saddle-point integration is done, the expression for the iVth trace reads 

jJ2 N R* 
-e *-n=i 



N 



(2tt 
dctB 



\2N 



NxN\ 



(4.23) 



To proceed, one needs to calculate the determinant of Bn x n, which is a tridiagonal non-symmetric matrix with 
additional non-zero elements in the upper right and lower left corners as well. The calculation of such a determinant 
can be found in the appendix |b], where we show that the dete r minan t is t he dif ference between the traces of two 
different products of 2 x 2 matrices. Combining equations ( 4.22| ), (B14) and ( 4.23| ) we are lead to 



trP- = £ 



p.p. \^ r Yil=N ^ 



(4.24) 



The inverted order of the indices at the product indicates that the matrices are ordered from left to right according to 
decreasing indices. The matrix My in the denominator is already the monodromy matrix for the purely dissipative 
part, M { d l) 



diag(^=p^-, 1). Unfortunately, a corresponding statement does not hold for Mi which is given by 



Mi 



1 



dui 







(4.25) 



Indeed, the monodromy matrix AfW for the entire Zth step (cf. (4.15)) has the form 



d Vl d m Si 



<9%Si) 



(du^SiY 



(dlSiKd^Si) -8%,Si 



2 



(4.26) 



when expressed in terms of derivatives of Si and /id- The fact that in M; both the indices I and I + 1 appear makes 
it even impossible to find a similarity transformation independent of I that transforms Mi to M^ l \ Nevertheless, we 
will show now that the traces of riz=iv % and rii=jv are ec l ua l- To do so, let us first pull out the common factor 
-1/dp^Si from M t and AfW by defining M ; = {-dufi^SfiMi and I7 (0 = {-dp^S^M^ . Then we separate the 



parts proportional to ^|=j from those independent of _ 



Mi 



dnd{vi) 
dui 



Ki+ Li, and M 



ji) = dn d (v l ) K m ■ L m 



dvi 



(4.27) 



We abbreviate ai = d^Si, j3i — d^Si and 7/ = dp l d IJll Si and find that the matrices on the rhs of the equations in 
( 4.27 ) have the structure 



K, = 



K (D = 



OLl + l 
1 

ai - 




Li = 



L« = 



01 




7f 



"7/ 






- OL101 0i 



(4.28) 
(4.29) 



Note that a/, /?/ and 7; are periodic in the index, i.e. ai + N = ai etc.. The product of the matrices Mi (or M''') gives 
rise to a sum of 2 N terms, each of which can be written as some factors multiplied with a trace over a product 
of matrices Ki and L/ (or K^> and L"\ respectively). To complete the proof we must show that the traces over 
products containing matrices Ki and L; remain unchanged if we replace all Ki by and all Li by . Note that 
these products are always ordered such that the indices decrease from left to right; and all indices from 1 to JV appear. 
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We will show the equality of the traces by establishing rules for calculating them for all possible combinations of 
A's and A's. The rules will turn out to be the same for both upper and lower index As and A's. 

Let us express our matrices in terms of spin-^ operators Si = o'j/2 and unity, where i = x,y,z and the <Ji are the 
Pauli matrices. In fact only raising and lowering operators S± = S x ± iS y , and ^ ± S z appear, 

Ki=S.+a l+1 (^+S z ) (4.30) 

L l = - 1 ?S++(3 l (±+S z ) (4.31) 

K® = -S+ + ai(^ + S z ) (4.32) 

L« = (j? - ai (3i)S- S z ) . (4.33) 

Such a decomposition greatly facilitates the calculation of products of A's and L's due to the simple multiplication 
rules of the spin operators. The rules can be summarized as 

S 2 _=S\= 0, S ± S T = (i ± S z ) (4.34) 

{\ ± S z ) n = (i ± S z ) for all natural n; S ± (± ± S z ) = (± ± S Z )S T = (4.35) 

(i±S.)5±=5 ± (iT5,)=5±. (4.36) 

Let us start by calculating the traces of products of lower-index matrices. 

1. Consider a product of matrices that contains only factors A;, Yli=j Ki where j > k are two arbitrary positive 
integers. With the above rules we easily find nf=j Kl = Kj Yli=j a i aim therefore, tr n!=jv ^ = IL^jv ai - 

2. A product containing only factors Li leads to Yli=j = 11?=/ 01 an d thus to tr n^=jv — Yil=N 

3. A product LfcAj gives L k Kj = (<x,-+i/3fc — lk)(k + ^ z )- With this result we can calculate a product which 
contains L's and A's, 



i+i i+i ! 

L k . . . LjKj-! .Ki = (J] AXajft - 7 /)( II + 5 *) > ( 4 ' 3? ) 



where k > j > i. Since (| ± S^)™ = (i ± Sy for all integer n, we can easily generalize the above equation to 
products containing several L and A blocks. We thus obtain the following rule for the trace of an arbitrary 
product starting with Ln and ending with K\\ 

• Replace all L/ by fti with the exception of the last Li in each L block. 

• Replace all A/ by a; with the exception of the last K\ in each A block. 

• Each interface between an L block and a A block (in this order, i.e. a term LjKj-\\) gives rise to a factor 
oijfy - 7?. 

We will call these three rules the standard rules. 



4. The next case to be considered is a prod uct st arting with Ln and ending with L\, with an arbitrary number 



of Li and K in between. We can apply (4.37) to all blocks with the exception of the last L block, to which 
we apply directly the rule for a pure L block. Besides the previous replacement of A; by ai and Li by (3i, one 
obtains a factor (i + S z )Li = L\. Since trAi = (3x the standard rules have to be modified only in as much as 
trijv • ■ ■ LK ... A' . . . Li gives rise to an additional factor Pi . 

5. Consider now a product starting with K^ and ending with L\. Due to the cyclic property of the trace we 
can bring the last L block to the left (e.g. trA^r ■ ■ ■ KjLj^i . . . L\ = trXj_i . . . L\Km . . . Kj), thus creating an 
interface L\Kn', and the product of matrices now starts with an L and terminates with a A matrix. Thus, the L 
and A blocks and their interfaces can be replaced according to the standard rules. Due to the cyclic properties 
of K\ and Li, also the interface L\Kn is of the type LjKj-i with j = 1 and is therefore to be replaced by 
oti(3\ — jf. In summary, the standard rules apply, but we get an additional factor — ~ff. 
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6. The last situation is a product starting with Kn and ending with K\. Here it is uselul to treat the first block 
Kjq . . . K j separately according to the rule for a pure K block and the subsequent sequence . . . K\ according 



to (4.37). One therefore encounters a factor tri^^(| + S z ) — cti. So again the standard rules apply, but we get 



an additional factor a\. 

Let us now consider the corresponding expressions for the matrices with upper indices resulting from the monodromy 
matrices. 

1 . A product containing only factors gives an expression with the same structure as a product of factors Li , 
Ul=j L[l) = L(k) U.l=j Pi- lt follows immediately that tr JILty = I\Ln A- 

2. A product consisting only of factors behaves differently from its counterpart with lower indices, for the 
last and not the first survives, 

Ui=j KW = K(k) I\t=j a i- However, since, trifW = ot\ we still have 

trnU^ (,) =nU««- 

3. A product of KM and is now simpler with on the left, K^L^ = (otifa - 7?)(§ + S z ) - p t S±. We can 
combine the rules for the pure K blocks, the pure L blocks and the ATL-interface to derive the product 



A (fe) ...A (j) L (j ~ 1) ...L w = 




(aiA-T?)(i + S.)-A^ 



(4.38) 



If we have several such KL blocks, each of them gives rise to a factor corresponding to the one in angular brackets 
on the rhs of the expression. But note that — J%){\ + S z ) — fliS + ] [{aifii — lf){\ + S z ) — /3/S+] = (aifli — 

7 2 ) \iaiPi — 7 2 )(^ + S z ) — &iS+ \ , i.e. from all KL blocks only the last angular bracket survives whereas all others 
are replaced by factors {pufii — 7 2 ) where the index is always the last index in the L blocks. We thus recover 
the rule of the replacement of LK interfaces by factors (aj/3j — 7, 2 )- In an arbitrary product starting with K^ N ' 
and ending in the last angular bracket leads to the factor (ai/3i — 7 2 ); and all K blocks, L blocks, and LK 
interfaces have to be replaced by the standard rules! Thus, we get just the same expression as for a product 
starting with Kn and ending with L%. 

4. In a product starting with K^ and ending with K^ it is useful to treat the last K block separately. If I is the 
last index of an L in the product, one is lead to a factor tr((aj/3; — "ff){\ + S Z ) — PiS + )K^ = tr(ai(3i — "/f)K < > 1 ) — 
(ai(3i — 7, )ai. Thus the standard rules apply, but we get an additional factor a±, in agreement with the result 
for the corresponding lower-index expression. 

5. In triW . . . L^K^-V ...KL... L^ we treat the first L block separately. One easily sees then that again the 
standard rules apply up to an additional factor j3\ , as was the case for the corresponding expression with lower 
indices. 

6. The last product to be considered is one that starts with £W and ends with K^\ We use the cyclic properties 



of the trace to shift the first L block to the right. Then ( 4.3S ) applies, and we see that the last L^> term on the 
right gives the same factor {ctjftj — t 2 ) as did the LjKj-i interface in the corresponding expression with lower 
indices. The standard rules apply without restriction. 

So we have shown that an arbitrary ordered product of factors K^> and L^ (I decreasing from A to 1 from left to 
right) has the same trace as the corresponding product in which all K^ 1 ' are replaced by Ki and all L^ by L\. This 

proves that trJTLw M i = ^Yl^N mW - 

The last thing to consider is the sign of each saddle point contribution for arbitrary A. We will show now that it is 
al ways p ossible to choose all minors of Q2NX2N real and positive. Therefore the sign of each saddle point contribution 
in ( 4.17| ) is positive. 



Observe first of all that without regularization all minors Di of Q2NX2N with the exception of the determinant 
dct Q2NX2N itself are zero. This is obvious for I — 1 ... N , since there the corresponding matrix is part of the upper 
left zero block of det Q2Nx2N- For I = N + m > N note that D[ contains a A x N upper left block which is zero, and a 
N x m (/z, rf) block in the upper right corner. Upon expanding Di after the first row one encounters subdeterminants 
with a (A — 1) x A upper left zero block and a (N — 1) x (m — 1) upper right (/i, rf) block. Both blocks together have 
A — 1 rows in each of which only the m — 1 elements on the right can be different from zero. Therefore the N — 1 first 
rows are always linearly dependent, unless m — N, the case that amounts to the full determinant. Thus, all minors 
Di with 1 < I < 2A - 1 are zero. 
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Suppose now that we add to a small quadratic term that vanishes and has a maximum at the saddle point 

0, but allows us to determine 



((j,, rf) = (/x sp ,0), i.e. a function — eX^iCCMi ~ Mspi) + Vi) with infinitesimal e > 0. If the original integral is 
convergent, the small addition will not change the value of the integral in the limit e 
the phase of all minors £>;. The matrices D; corresponding to Di get all replaced by DJ = D/ + el; 



unit matrix in I dimensions. For 1 < I < N we are immediately lead to Di 
expand D[ in powers of e and get D[ = D t + etrD; + 0(e 2 ) = etrD; + 0(e 2 ) 
the second derivatives in the (ry, rf) block of Q2Nx2N- They are given by 



where 1; is the 
= e' > 0. For N + 1 < I < 2N - 1 we 
To determine the traces trD; we need 



d m d m ^N = 



d Vk R(nk,Vk+v,r)k) 

32 



{(h k d^ k S k f 



2 



+4 — ^ — — ^ — Ok-i,l + 4 a2 u/ — -Ok+l,l 



d v k R ^k-i,Vk]rik-i, 



(4.39) 



Since 9^ fc i?(/ifc, Vk+i\ Vk) < and d^ k R(fik~i, ^k', Vk~i) < for all k at the saddle point, the diagonal elements of the 
(r], rf) block of Q2NX2N are all real and positive definite (remember that we pulled out a minus sign in the definition of 
Q2N x2N in terms of second derivatives). Therefore trD; is real and larger than zero for all I. Since also det Q2./Vx2jv 
is real and positive, t he sig n of each saddle-point contribution is now determined to be positive (cf.(A3)). The proof 
of our trace formula (4.17) for the TVth trace i s her ewith complete. It is worth mentioning that the comments made 
abo ut th e decoherence built in in the formula ( 3.11 ) for trP, the limit r — > 0, and the generality of the formula apply 
to (p7|) for trP^as well . 



V. EVALUATION FOR AN INTEGRABLE DISSIPATIVE KICKED TOP 

We now propose to evaluate explicitly trP N for our integrable dissipative kicked top. To do so, we need the periodic 
points of the classical map P c \ and their stabilities. 



A. Periodic points 

Let us start by determining the fixed points of P c i- It is most convenient to work with Cartesian coordinates, 
x = J X J J, y = Jy/ J, and z — J z / J — /1 (with x 2 + y 2 + z 2 = 1). The coordinates x' N , y' N and z' N immediately after 
the iVth rotation are related to those immediately before the iVth rotation (xn, yN and z^) by 

x'n = X N cos/3 + zjq sin/3 (5-1) 
Vn = Vn (5.2) 
z' N = — xn sin/3 + zjy cos/3 . (5-3) 

The dissipation takes these intermediate coordinates to those immediately before the (N + l)st rotation, 

W+. = . ■ . (5-5) 

cosn t — z N smn r 

z N+1 = - s[lihT + z N coshr _ ( 5 6 ) 

cosh t — z' N sinh r 

For r — 0, the dissipative part becomes the identical map and the fixed points are those of the rotation, thus 
xat+i = xat = (0, ±1,0) for all (3. 

For t ^ 0, the equation for the y component demands yj\r — or coshr — z' N sinhr = 1. Let us consider the latter 
equation first. It leads to the fixed points 



[3 t 1- cos/3 -2tanh 2 § r . 

x N = -cot-tanh-, y N = ±\ , ZAr = -tanh-. (5.7) 

2 2 y 1 — cos p 2 



15 



They can exist with real um only if sinh 5 < | tan 5|. If we compare (5.7) with x t+ i = = (0, ±1, 0) at r = 0, we see 



that the pair of fixed points shifts downwards in z and towards the y = plane as r increases. When t — t c given by 
sinh y = | tan the two fixed points coincide in the plane y — 0; a bifurcation takes place and one has to consider 

yN = 0. 

This case is treated most easily in polar coordinates, since yN = means = and x depends only on the polar 
angle 0. The rotation becomes trivial, 0' = On + 0, w here should be allowed to run through — tt...tt in order to 
avoid <p ~ tt. The overdamped-pendulum trajectory ( |2.9| ) implies 

tan " N+1 — e T tan — = e T tan - — — (5.8) 
2 2 2 y ' 

throughout the interval = —n...n. With the substitution un — tan ^M- the periodic-point condition can be written 
as a quadratic equation, 

tan- u 2 N + u N (l - e T ) + tan^e T = 0. (5.9) 

The two solutions 



u Nl . 2 = I - sinh - ± ^sinh^ - - tan 2 '- ] (5.10) 



T , /_._,_a r . o 



are real iff sinh ^ > | tan ^|, i.e.r > r c . Transforming itjv back to /x, we get 



u = ± cos —\ 1 — (coth — sin — ) 2 — coth — sin 2 — . (5.11) 
' 2 V 2 2 y 2 2 y J 

For r — > oo, the lower one of the two fixed points moves to the south pole (J, = — 1 for all values of [3. The upper one 
moves to /i — — cos/3, such that the inte rmed iate point immediately after the rotation is at the north pole, ji = 1. 



The final momentum v is obtained from (5.11) by reversing the sign in front of the last coth-| 



The boundary between classically allowed and forbidden rotation is given by 0^r + 0' = (3 (see section [I A ), or 

sin 2 /?- y 2 -v 1 + 2yu> cos = 0. (5.12) 



One easily checks that the two fixed points given by (5.11) have the remarkable property that they lie exactly on this 



boundary, f or a ll values of t > t c and 0, as long as the fixed points exist. Unfortunately, in this case the semiclassical 
expression ( |2.4| ) for the unitary propagator ceases to be valid. To avoid this problem we will restrict ourselves to 

T < T c . 

We show in appendix |c| that has the same and only the same fixed points as P c i . Therefore the restriction r < t c 
is sufficient as well for all higher traces. We also remark that N0 should not be an integer multiple of 27r, otherwise 
already the semiclassical approximation of the unitary part breaks down ]l£f| . 

In FIG|l] we show phase-space portraits. The northern and southern hemisphere of the spherical phase space 
spanned by the variables /i and <f> were projected on the equatorial plane /i = 0. The two symmetry related fixed 
points are clearly visible for r < t c as elliptic fixed points, whereas for r > r c they become a pair of a point attractor 
and a point repeller. 



B. Stability matrices 

The only non-trivial element of the monodromy matrix Md for the dissipative motion is given by — = 

\Z„'i , whereas d ^ d g^'^ — 1 and the two other elements are zero. For the fixed point at t < r c , v — —pi, even this 
non-trivial element equals one, and thus Md — 1. In the vicinity of its two fixed points our dissipative map is actually 
area-preserving . 

The elements of M r can be found from differentiating the equations of motion for the rotation written in /z, <f> and 

y! = v, ft, 

\J\ — v 2 sin <fi' = \/\ — sin <j) (5.13) 
v = fi cos — \]\ — cos cf) sin . (5-14) 



1G 



We then obtain at v = —\i 

ih ' sin 2 /? . , 

; sm (p (5.15) 



1 + cos j3 
sin 2 

.^/3+— — 
O/i 1 + cos p 



df 1 ' i , n sin 2 /3 2 , 

cos/3+- --rcos 2 0, (5.16) 



and thus at the fixed points ([3.7]), 



trM = trMi = „— (tanh 2 - + cos ,5) . (5.17) 

1 — tanh 5 2 



If we insert this and trM^ = 2 in (3.11), we obtain the final result for r < r c , 

1 - tanh 2 ^ 

trP = 2 -^- . 5.18 

1 - cos 0- 2 tanh 2 \ 

A factor 2 has canceled since two orbits contribute with equal weight. Clearly, trP diverges at the bifurcation. 
Since M\ is unitary for ^jf- — 1, its eigenvalues can be written as A1/2 = e ±lS . We have therefore trMi = 2 cos 5 and 

trM^ = cos(NS) = cos(7V arccos trJ ^ 1 ), or 

1 i i\t /tanh 2 f-+cos/3ss ' (5.19) 

1 - cos(iVarccos( z )) 



C. Comparison with numerics 



In FIG.|| we have plotted the first trace as semiclassically approximated by ( 5.18| ) together with the numerically 



obtained quantum counterpart for different values of j. The agreement is good for r > 1/J and r not too close to the 
critical r c , where the bifurcation takes place and the semiclassics diverges; the agreement improves with increasing 
j. The rather erratic behavior of trP as a function of J when r < 1/J is a signature of the importance of quantum 
effects in the dissipative part. 

In FIG.|| we compare the semiclassical and quantum mechanical (j = 20) result for tiP N as a function of N. The 
agreement is rather good for the first about 10 traces. Both results show oscillations, and the semiclassical result 
reproduces perfectly their period and phase over the entire N range. However, the exact traces decay exponentially as 

trP N N -^3° l ; whereas the semiclassical approximation ( 5.1 9| ) shows undamped oscillations. Note that all eigenvalues 



Aj of P are bound by |Aj| < 1 due to stability reasons. In general one has even |Aj| < 1 with the exception of one 

eigenvalue Ai — 1 which corresponds to a stationary state. Since trP w = X)i=i +1 '' ^ the exact traces must indeed 
decay exponentially towards 1. 

We have studied the decay of the exact traces numerically in more detail. It is an exponential decay of the form 

trP-W ~ 1 + const.e J , where /(r, (3) is some function of the parameters. Thus, the effect of the dissipation 
is very small and vanishes for fixed finite N in the classical limit. Its onset should be visible one order beyond our 
present semiclassical approximation. Note that the absence of classical periodic orbits with phase-space contraction 
alone is not enough to expect a spectrum corresponding to a unitary map (i.e. eigenvalues on the unit circle). This 
is due to the fact that beyond the classical information entering via the periodic points our trace formula always 
contains a dissipative effect on the quantum mechanical level, namely the complete destruction of coherence manifest 
in the single sum over periodic points and the cancellation of the actions from the unitary part. 

In order to arrive at the precision at which the exponential decay should be seen, it would be necessary to obtain at 
least the purely dissipative propagator D to accuracy 0(1/ J). So far our semiclassical calculations of trP have been 
restricted by the requirement 1/ J <C t <C J which marked the limit of validity of the semiclassical approximation of 
D. For the calculation of trP^ one should have correspondingly N/J <C Nt <C J. 

In order to reconstruct the whole spectrum one would need traces from N = 1 to N = (2J) 2 , such that r could 
only be less than 1/J to fulfill the last condition. But this conflicts with r ^> 1/J. We conclude therefore that 
without pushing semiclassics one order further, it will never be possible to calculate a whole spectrum. However, 
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this might be a too ambitious goal anyway: Even for moderate dissipation (r ~ 2 for example), most eigenvalues are 
exponentially small and can not be reconstructed numerically even if one knew the exact traces! We tested this by 
calculating traces from exact quantum mechanical eigenvalues. From the traces one can get the coefficients of the 
characteristic polynomial via Newton's formulae and therefore in principle back the eigenvalues. But the polynomial 
becomes exponentially small for its argument close to zero, which makes it virtually impossible to accurately determine 
its roots. For large values of J one faces moreover the problem that the higher traces (N ~ J 2 ) are all equal to unity 
up to tiny corrections of order C(exp(— J)) and therefore carry no more information. 

On the other hand, it is quite likely that in typical experimental situations only the slowest decaying eigenmodes and 
corresponding eigenvalues can be measured. Future investigations will have to show, in how far these can be obtained 
from a limited number of traces. 



VI. CONCLUSION 

We have established a semiclassical periodic-orbit theory for a dissipative quantum map with an area-nonpreserving 
classical limit. We have shown that traces of arbitrary integer powers of the propagator of the density matrix can be 
written as a sum over classical periodic points. Our results generalize the well known formula due to Tabor for area- 
preserving maps and are in the spirit of Gutzwiller's trace formula. Both the exponent and preexponential factor in 
Tabor's formula are modified. A "diagonal approximation" which amounts to neglecting interferences between 

different orbits in suitably averaged quantities arises automatically for our dissipative case: Interferences between 
different orbits are destroyed by the dissipation. Therefore the exponent in the trace formula is purely real. Only the 
action from the dissipative part of the map remains. Even though only classical information is needed to evaluate our 
trace formula, the quantum mechanical effect of the destruction of coherence by dissipation is automatically built in. 
In the preexponent the total monodromy matrix, but also the monodromy matrix of the dissipative parts alone enter. 
The Maslov index is always zero. 

Our theory allows to quantify the effects of the dissipation on the spectrum of the propagator even when dissipation 
is so strong that it not only leads to dephasing but also to a modification of periodic orbits. Our example of a 
dissipative integrable kicked top showed, however, that even though the classical map may be area-nonpreserving 
globally, it might have periodic orbits close to which it is area-preserving. If only such orbits exist, the only influence 
of the dissipation on the spectrum is via the destruction of interferences between different periodic orbits. 
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APPENDIX A: SADDLE-POINT METHOD FOR A COMPLEX FUNCTION OF SEVERAL VARIABLES 

Let x represent M real variables x\ . . . xm, and let F(x) and G(x) be complex valued functions of x with 5RF < 
in the volume V in 1Z M over which we are going to integrate. Suppose that V contains a single non-degenerate 
stationary point Xo with d Xi F(x.o) = 0, i = 1 . . .M. Let us denote by Qmxm the matrix of the negative second 
derivatives, (QMxM)ik — —d Xi d Xk F(x) taken at x = xo. The condition of non-degeneracy means det Qmxm =/= 0. 
We then have for J — > oo 

[ d M xe JF ^G(x) = G(x )J .., ^™ e Jf(xo) -* IndQ ™(l + 0(l/J)). (Al) 

Jv V J \detQ M xM\ 

Here IndQjvfxM is the index of the complex quadratic form x of M real variables, 

M 

X = X! (QMxM)ijXiXj ■ (A2) 

The index is defined via the minors Dk = det ||(Q.wxM)y||) 1 < i, J < As of Qmxm as 
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IndQ 



MxM 



^2&rgPk, - 7T < argp fc < vr. 



fe=i 



pi = D\ = (Qmxm)u, Pk 



D k -x' 



k = 2...M. 



(A3) 
(A4) 



The restriction — 7r < argp^, < ir on the phases fixes uniquely the overall phase of the saddle point contribution. 
Without this restriction IndQi/xAf would only be defined up to multiples of 2-7T, and this would lead to an overall 
phase ambiguity corresponding to the choice of sign for the square root in (Al). 

In the case of a Dk = we add a term —e(xi — ^oj 2 to the function F(x). Such an addition does not change the 
convergence of the integral (in the limit of J — *• oo), nor its value for e — > 0. However, such a small term may bring 
the Dk awa y fro m zero and therefore allow us to determine its phase. 

Formula (Al) has the structure familiar from the SPA of a complex function of a single variable. The term 
— jlndQjwxM leads to a phase analogous to Maslov's ±7r/4, but that phase can now take on any value between 



APPENDIX B: THE DETERMINANT OF A TRIDIAGONAL, PERIODICALLY CONTINUED MATRIX 



The N x N matrix whose determinant we want to calculate has the structure 

/ o u mi ... 

0-21 022 a 2 3 . . . 

o '■. '■. '•■ 



A = 



am \ 




\ 1JV1 







(Bl) 



0>N,N-1 a NN J 



We first consider the simpler problem where a±N = ajvi = 0; the corresponding matrix will be called Aq. The 
determinant of Ao can be expanded along the diagonal by using the block decomposition rule for determinants |p8| . 
If A is decomposed into smaller matrices W (n x n), X (n x to), Y (to x n), and Z (to x to, where n + in = N) and 
if the inverse of W exists, we have 



det 



W 
Y 



detWdet(Z -YW^X) 



(B2) 



We apply the rule to det Aq by choosing W — an. Then Z = A^ is the matrix in which the first column and row 
are missing. In the matrix YW~ X X only one element is different from zero, namely (YW^ 1 X)ik = t>i,2<5fc.2 a21 ° 12 . 
Note that we keep the numbering of indices of the original matrix A. Thus, in the determinant on the rhs of |l32| ), 
which is now only (N — 1) x [N — 1), the indices run from 2 to N . The new matrix whose determinant has now to 
be calculated differs from the corresponding block in A just in the upper left element d22, and is tridiagonal again. 
The procedure may therefore be iterated. In each step the dimension of the matrix is decreased by one; the only 
matrix element changed is always the one in the upper left corner. It gets replaced by the original one in Aq minus 
the product of the two secondary diagonal elements flanking it divided by the upper left element of one step earlier. 
One obtains therefore a product of renormalized diagonal elements, 



N 



det Aq = TT cj 

3 = 1 



with ci = an and the iteration law 



a j,j-i a j-iJ 

Cj-l 



(B3) 



(B4) 



If we evaluate the product starting from the last term, we find that all Cj's with the exception of the one with the 
smallest index cancel, such that 



N 



9n 



(B5) 
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The functions /„ and g n depend only on the matrix elements, not on the c/s. In the notation of (B5) det Aq is given 
by det = ci/i + <7i, and the final condition for /„ and g n reads /jv = 1, 9n = 0. A recursion relation for /„ and g n 



is easily found by equating rii = n-i 



l(Cnfn + 9n) = C„_l(a n „/ n + Q n ) 



We obtain a linear relation which we write with the help of the matrix 



Qr 



Q"nn 1 
Q>n,n— l^n— 1,ti 



-l&n— l,n fn — ^n—xfn—X ~i~ Qn- 



(B6) 



' fn - 1 )=Qn( fn 

\9n-iJ \g n 

The linear combination C1/1 + g\ at our interest therefore reads 



(B7) 



(ci,l) 



(an, 1) 




The product can be extended down to j — 1 by noting that the vector (an, 1) is just the first row of the matrix in the 
product that would carry the index 1. In the second row the product of elements —121,000,1 appears. Its value is at 
this stage completely arbitrary. We define it by periodically continuing the indices, i.e. a^o = di,N an d <Xo,i = <2jv,i- 
We then arrive at the result for the determinant of the tridiagonal matrix Aq, 



det A, = (1,0) 




(B8) 



Consider now the case where ffli,jvttiv,i 7^ 0. The problem of calculating det A can be reduced to the problem of 
determinants of purely tri-diagonal matrices by expanding det A after the first row and the resulting submatrices after 
the first column. One eliminates in this way the rows and columns with the disturbing additional corner elements. A 
matrix appears, which is defined as the matrix A from which both the first and Ath rows and columns have 

been taken out. The reader will easily verify that 



det A = det A Q ~ a 1N a N1 det A (1,JV ) + (-1 



,iV+i 




(B9) 



For clincini = the expression reduces to det Aq as it should. It will be useful to rewrite the sum of the two products 
as trace of a diagonal 2x2 matrix which is the product of matrices Rj = diag^j-i, cij-ij). Let us apply now our 
result for tri-diagonal determinants (BE). We obtain 



det A = (1,0) 




CL\,N a N,l 




+ (-1) 



N+l 



N 

txII R 



Now observe that 



(0,l)Qi 




Q 



iV 



cti,iv<2Ar,i(l, 0) and 
0j 



(B10) 

(BH) 
(B12) 



Thus, while the first term in (B10) gives the upper left matrix element of the product of matrices, the second one 
gives the lower right. We therefore have 



/V 



detA = trJjQj + (-l) 



N+l 



N 

trf]^ 



(B13) 



The order of the products is such that they start with the matrices with index 1 on the left. Given the ordering of 
the monodromy matrices that appears in the main text it is useful to inverse the order. This can easily be done by 
taking the traces of the transposed product, which leads to the final result 
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det A = tr J| 

j=N 



"JO 

1 







+ (-1) 



N+l 



trn 

j=JV 



dj-ij 



(B14) 



Do not apply the formula for N = 2, since we assumed additional elem ents in the corners. For N = 2 the corners 
collapse with the secondary diagonal and are therefore over counted in ( B14 ). The lowest dimension for which the 
formula works is N = 3. 



APPENDIX C: FIXED POINTS OF Pg L 



Let us first consider the case y n ^ 0. Equating un+i with yi in (|5.5|), one is lead immediately to n^li( cosrir — 
sinhrZj') = 1. Since the same denominator also appears in Xn+i and Zn+i when expressed as functions of x\ and z\, 
the remaining map (x\, z\) — > (xjv+i, 2zv+i) is in fact linear! Foreseeing this we can build it up by iterating a linear 
map A which acts on the two dimensional vector = (xi,zi) as 



0+i = AO + b . 
The matrix A and the vector b are explicitly given by 



A = 



cos/3 



sin/3 



cosh t sin /3 cosh r cos (3 



b = 





— sinh t 



(CI) 



(C2) 



The condition for the fixed points of P c i reads (1 — A)0 = b. The determinant of (1 — A) is (1 — cos/3)(l + coshr) 
and therefore non-zero as long as there is any rotation at all. The equation is solved easily and gives back the fixed 
point at Z\ = — tanh but additionally we now see directly that it is unique in Xi, Z\, as the map is linear and the 
determinant non-zero. The y component is obtained from normalization, yi — ±yl — x\ — z\. These are the two 
basic symmetry related fixed points of P c i. 

The fixed points of P\ follow from = 0j = A( 2 +b = A 2 + (1 + A)b, i.e. Ci = (1 - A 2 ) _1 (l + A)b = (l-A)" 1 ^ 
One easily verifies by complete induction that the formula for the iVth iteration is 



TV 



(l-A«)( 1 = £A«)b^ 



I- A 



N 



-A 



(C3) 



and therefore again just (1 — A)0 = b, the equation for the fixed point of P c ; 



leads exactly to the same 



For (j) = the situation is similar: One easily verifies that un+2 — un m equation (5 
equation ( |S.9[ ) as did un+i = ujsr, and the same is true for all higher iterations. This completes our prove that P c 
has the same and only the same fixed points as P c i. 
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FIG. 1. Phase-space portraits of the classical map of the dissipative integrable kicked top [fi = 1.2). We projected the 
northern and the southern hemisphere (n > and /i < 0, respectively) of the spherical classical phase space on the equatorial 
plane fi = 0; (a) 0.6 = r < r c ~ 1.28, the two fixed points are elliptic fixed points; (b) 1.3 = r > t c ~ 1.28, the fixed points 
form a point attractor/repeller pair. The repeller in the northern hemisphere can be made visible by iterating the inverse map. 



FIG. 2. 

Comparison of semiclassical (dash ed lin e) and quantum mechanical result for trP as a function of the damping r 
at (3 = 2.0. The semiclassical result (|5.18 ) diverges at r = r c ~ 2.45, and the quantum mechanical trace follows this 
divergence the further the higher the value of J (circles j — 10, squares j — 20, diamonds j = 40). For r < 1/J the 
problem is purely quantum mechanical and our semiclassical treatment seizes to be valid. 



FIG. 3. 

Comparison of semiclassical and quantum mechanical result for trP^ ((3 = 0.5, r = 0.2) as a function of N. 
FIG. (a) for j — 20 shows that the agreement is fairly good for about the first 10 traces and the quantum mechanical 
oscillations are well r eproduced over the entire N range. However, the exact traces decay, whereas the semiclassical 
approximation ( 5.19 ) only shows the oscillations (same symbols as in FIG.l). In (b) the exponential decay of the 
exact traces is studied numerically for different values of j. It is of the form trP N ~ 1 + const. exp(^|j) for (3 = 0.5 
and t = 0.2. 
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